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SUMMARY 

A new computational simulation method is presented to evaluate the dynamic 
aspects of composite structural response and durability that have not been sim- 
ulated previously. Composite structural behavior under any loading condition, 
geometry, composite system, laminate configuration, and boundary conditions can 
now be simulated. Structural degradation, delamination, fracture, and damage 
propagation are included in the simulation. An angle-plied composite plate 
structure under normal impact loading is used as an example to demonstrate the 
versatility of the simulation method. 


INTRODUCTION 

The behavior of fiber composite laminates during progressive damage and 
fracture has become of increasing interest in recent years due to the multitude 
of benefits that composites offer in practical engineering applications such as 
lightweight airframes, engine structures, space structures, marine and other 
transportation structures, high-precision machinery, and structural members in 
robotic manipulators. Composite structures lend themselves to tailoring to 
achieve desirable characteristics such as a high strength to weight ratio, 
dimensional stability under extreme thermal and hygral fluctuations, and the 
capability to allow controlled detectability such as in the Stealth technology. 
Because of the infinite possibilities with material combinations, composite 
geometry, ply orientations, and loading conditions, it is essential to have a 
reliable computational capability to predict the behavior of composites under 
any loading condition, geometry, composite material combinations, and boundary 
conditions. The predictions of damage initiation, progressive damage, and pro- 
gressive fracture are important in evaluating the safety and reliability of 
composite structures. 

The behavior of composites during progressive fracture has been investi- 
gated both experimentally and by computational simulation (refs. 1 and 2). The 
concept and computational foundations of CODSTRAN (Composite Durability STRuc- 
tural ANalysis) were first laid out by Chamis and Smith in 1978 (ref. 1). The 
computational simulation of progressive fracture predicted by CODSTRAN was 
found to be in reasonable agreement with experimental data from coupon tension 
tests conducted at NASA Lewis Research Center (ref. 2). From the very begin- 
ning, CODSTRAN was designed to be a modular and open ended computer code. The 
modular nature of CODSTRAN has enabled the upgrading of nearly the entire code 
as new and improved modules have become available for composite mechanics, 
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cTnV»/ Ura * anal y s ! s .' an( * lure criteria. Due to the open-endedness of COD- 
STRAN, recent additions to the code have enabled monitoring the variations in 
structural properties such as natural frequencies, vibration mode shapes, and 
buckling modes during progressive fracture (ref. 3). In these past studies 
computational simulation of progressive fracture and failure has been primarily 
concentrated on the effects of monotonical ly increasing static loading on the 
progression of damage and fracture. The objective of the current research is 
to develop and demonstrate the computer simulation of composite behavior with 
damage progression and fracture under dynamic loading applied normal to the 
laminate, including loading rates that are sufficiently abrupt to be considered 
impulsive . 

For the purpose of this discussion, composite structural damage under 
dynamic loading may be classified into three levels: (1) Initial Damage, (2) 

Damage Growth, and (3) Fracture Progression, (l) Initial Damage refers to the 
initiation of structural deterioration. It is usually localized under the 
applied normal dynamic loading and may involve one or more plies of the lami- 
nate. This is generally a matrix dominated failure (matrix microcracking 
caused by transverse as well as shear stresses within the damaged plies). (2) 
Damage Growth refers to additional damage at the initial damage locations with 
increased loading. There may be additional failure modes and/or additional 
plies participating in the damage at this secondary level. (3) Fracture Pro- 
gression refers to the growth of damage and fracture from a localized area 
under the impact loading to other parts of the composite structure. 


BACKGROUND 

CODSTRAN (Composite Durability STRuctural ANalyzer) is a computer code 
developed for the analysis of composite structures; keeping track of progres- 
sive damage resulting from the application of loads (refs. 1 and 3). It con- 
sists of an executive module, an analysis package for micromechanics and 
macromechanics (ICAN), and a finite element analysis module (MHOST) , The ICAN 
(Integrated Composites ANalyzer) (refs. 4 and 5) and MHOST (ref. 6) computer 
codes are originally independent codes that have been integrated into CODSTRAN 
as analysis modules. 

The ICAN computer code incorporates the cumulative result of composites 
research at NASA Lewis Research Center on polymer matrix multilayer angle-plied 
composites. ICAN is capable of determining the ability of the composite lami- 
nate to endure stresses and deformations caused by applied loading and environ- 
mental effects such as temperature and humidity. ICAN predicts all possible 
internal damage and fracture in the composite laminate. Fourteen distinct 
failure modes are checked by ICAN at each load increment during a CODSTRAN 
analysis stage. These failure modes include the failure criteria associated 
with the positive and negative limits of the six stress components, a modified 
distortion energy fai lure cri terion , and delamination due to relative rotation. 
Another important capability of ICAN is to compute the anisotropic elastic con- 
stants from the physical information on the composite laminate. ICAN utilizes 
a resident data bank that contains the properties of typical fiber and matrix 
constituent materials , wi th provisions to add new constituents as they become 
available. The piecewise linear constitutive relationships determined by ICAN 
are supplied to the MHOST finite element analysis module for the definition of 
nodal finite element properties at each load increment. The executive module 
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of CODSTRAN determines the appropriate step size for each load increment and 
directs the computational simulation of nonlinear-inelastic, large deformation 
progressive damage and fracture behavior under loading. 

The integration of ICAN with a general purpose structural analysis code 
results in a piecewise linear incremental loading analyzer with equilibrium 
checks at each local iteration for each load increment. The current version of 
CODSTRAN uses MHOST as the structural analysis module due to the ability of the 
MHOST quadrilateral shell element to accept the composite laminate force- 
deformation relations predicted by ICAN. 

As an example of the capabilities of the CODSTRAN code existing prior to 
the current research, the load-displacement relationship computed by CODSTRAN 
for a typical composite plate structure is presented in figure 1. This struc- 
ture is loaded by a uniform static tensile loading applied in the longitudinal 
direction in the plane of the composite. The apparent softening of the simu- 
lated load-displacement relationship with increased loading is because of the 
initiation and progression of damage and fracture in the composite plate. 

After the application of the last load increment the plate is torn into two 
pieces . 

The CODSTRAN incremental loading procedure uses an accuracy criterion 
based on the allowable maximum number of damaged nodes during the application 
of a load increment. If too many nodes are damaged or fractured during a load 
increment, incremental loads are reduced and the analysis is repeated. Other- 
wise, if there is an acceptable amount of incremental damage, the load increment 
is kept constant but the constitutive properties and the structural geometry are 
updated to account for the damage and deformations in the previous increment. 

The structure is then reanalyzed for further damage and deformation. If, after 
an incremental analysis step, there is no damage as determined by ICAN, the 
structure is considered to be in equilibrium and an additional load increment 
is applied. Analysis is stopped when global structural fracture is imminent. 

During the computational simulation of progressive fracture, CODSTRAN also 
keeps track of degradation of structural integrity such as natural frequencies, 
vibration modes, buckling loads and buckling modes of the structure. CODSTRAN 
has been experimentally verified for the simulation of progressive fracture and 
failure under monotonical ly increasing static loading applied in the plane of 
the composite (refs. 1 and 2). 

In the present paper, CODSTRAN is augmented to conduct the simulation of 
progressive fracture under dynamic loading applied normal to the laminate. 
Loading rates for the newly augmented CODSTRAN code may vary from relatively 
slow loads that could be simulated by a pseudo-static analysis, to much higher 
loading rates that may be classified as impulsive. All dynamic effects are 
taken into account within the limitations of the thick shell finite element 
model. Effects that are not yet included are dynamic wave propagation across 
the thickness of the composite laminate and strain rate effects in the consti- 
tutive relationships. The dynamic version of CODSTRAN uses a similar incremen- 
tal loading procedure as in the static version. However, inertial effects are 
also taken into account. The dynamic time history is simulated by using the 
Newmark implicit integration algorithm that is available in the MHOST module. 

It is necessary to maintain the incremental loading process while updating the 
nodal coordinates, velocities, accelerations, and constitutive relations in the 
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dynamic analysis procedure in order to closely simulate the realistic progres- 
sion of damage and fracture. In the augmented CODSTRAN program the effective 
loads at each dynamic load increment are expressed in terms of the incremental 
loads and the structural mass, nodal accelerations, velocities, and 
displacements. 

The dynamic loading may be specified as any arbitrary function of time. 

The time domain is subdivided into a sufficient number of time intervals for 
accurate dynamic analysis. High frequency vibration modes are usually signifi- 
cant during dynamic impact loading, requiring a relatively small time step 
size. However, structural geometry and laminate properties do not change sig- 
nificantly at each time step. The number of time steps to cause considerable 
c ange in structural geometry and/or to cause possible laminate damage consti- 
tute the dynamic load increment. The load increment size is dictated by the 
criterion on the allowable number of nodes of the structure with incremental 
damage. The duration of a load increment usually takes many time increments. 
Laminate properties and structural geometry are updated at the end of each load 
increment as computed. Within a load increment the equations of motion are 
solved to obtain the time history response using the number of time increments 
contained in that load increment. Time histories for consecutive load incre- 
ments are chained together with the appropriate continuity conditions to obtain 
the overall time history response. The temporal discretization requirements 
depend on the rate of the dynamic loading application as well as on the struc- 
tural properties. The spatial discretization needs are mainly dependent upon 
the required extent and accuracy of damage and progressive fracture analysis 
after the initiation of damage and large structural deformations. 


EXAMPLE COMPOSITE STRUCTURE AND ANALYSIS 


As an example of the augmented simulation capabilities in CODSTRAN an 
intermediate stiffness T-300/Epoxy (Thornel-300 graphite fibers in an epoxy 
matrix) composite structure is selected for analysis of dynamic response under 
impulsive loading A composite square plate (3- by 3-in.), supported along 
its tour edges is loaded by a concentrated normal dynamic load at the center 
tu 1 ®*- 2 ' ‘ ,8 e deformations are included in the computational simulation. 

The finite element model consists of 64 plate elements as shown in figure 3. 

The finite element model is deliberately chosen to be somewhat coarse to ensure 
a reasonable computer turnaround time for this demonstration example. The 1am- 
inateconf juration is +45/-46/-46/+4S/, with a total composite plate thickness 
or 0.02 in. All edges have been restrained against displacement in the z 
direction in all cases. Within this constraint, two types of support condi- 
tions are considered In the first type, support nodes are allowed to move in 
the x y plane (the plane of the undeformed composite plate) In the second 
type of support, boundary nodes are restrained against displacement in all 
directions. The first and second types of support conditions will be refer- 
enced to as released and restrained supports, respectively. 


The concentrated transient load is assumed to increase linearly with time 
as shown in figure 4, until global fracture of the composite plate structure. 
Ihree loading rates are considered. These loading rates are 0.1, 2.0, and 
3.5 lb/psec. Figure 5 shows the dynamic load-deflection histories at these 
three rates for the first type (released) boundary support condition. The dis- 
placements are very small for the first load increment because of the large 
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inertial load imparted when the structure is initially set in motion. Also, 
the higher the loading rate, the smaller the displacement because of the higher 
inertial force. After the first load increment, there is a dramatic increase 
in displacement, especially for the 0.1 lb/psec loading rate. In this case, 
the time is sufficiently long for the entire plate to be set in motion. After 
a center deflection of approximately 0.3 in. membrane forces become significant 
and the incremental deflections are reduced. It may be noted that even though 
the support nodes are free to move in the x-y plane, membrane effects are 
developed due to the restraint of boundary nodes in the z direction and the 
square geometry of the plate. If the support nodes were not restrained in the 
z direction, the corners of the plate would lift up and fold along the diago- 
nals to allow for a truly unrestrained deflection of the plate without membrane 
forces. The load corresponding to 0.3 in. center deflection is 60 lb. This 
loading level also causes initial local damage in the composite laminate 
(ply 1 - the first ply on the tension side of the plate). The initial damage 
is caused by the ply transverse tensile stresses reaching their limits. 

Figure 6 shows the contours of ply 1 transverse stresses that cause the 
damage. Figure 7 shows the much lower intensities of the same stresses after 
CODSTRAN accounts for the damage, redefines the material properties, and 
reaches a dynamic equilibrium with the redefined properties. Numerical results 
show that the transverse ply stress is zero at the damaged center node. How- 
ever, this zero ply stress does not appear as a contour in figure 7 because 
the damage is highly localized involving only the center node under the load. 
Figure 8 shows the deformed finite element mesh with displacement contours 
immediately after initial damage. The deformed structure continues to carry 
additional load, sustaining approximately constant incremental displacements 
with each load increment. Other plies participate in the laminate damage 
growth at approximately 240 lb loading. The maximum deflection at the center 
node under 240 lb loading is 0.7 in. It may be noted from figure 5 that there 
is no apparent change in the overall stiffness due to localized transverse 
damage in the plies of the laminate. This lack of sensitivity to damage is 
because membrane forces are much more significant compared to flexural forces 
under large deformations. Localized damage in the ply transverse directions 
reduces only the flexural stiffness of the laminate. The entire structure par- 
ticipates in the dynamic response similar to a static loading. The deformed 
finite element mesh with displacement contours corresponding to 240 lb loading 
at 0.1 lb/psec is shown in figure 9. 

For the 2.0 lb/psec loading case, damage is initiated by the same trans- 
verse stress limit criterion in the tension plies. However, initial damage 
occurs at a higher dynamic load of 75 lb. The center node displacement corre- 
sponding to the 75 lb damage initiation load is approximately 0.07 in., as it 
may be observed from figure 5. For the 3.5 lb/psec loading rate initial dam- 
age occurs at a load of 105 lb. The maximum deflection corresponding to 105 lb 
for this case is also approximately 0.07 in. The 2.0 and 3.5 lb/psec loading 
rates both cause the same type of response that is characterized by significant 
local deformations that do not engage the entire plate. 

Figure 10 shows the deformed finite element mesh and the displacement con- 
tours corresponding to the initial damage load of 105 lb applied at the rate of 
3.5 lb/psec. The initial damage mode is the transverse stress limit as in the 
slower loading cases. However, in this case three plies are failed simultane- 
ously. Figures 11 and 12 show the ply 1 transverse stresses, before and after 
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dynamic equilibrium, respectively, at the time of initial damage. Further 
loading after initial damage does not significantly change the displacement 
pattern. However, the local deformations under the load are further increased 
and longitudinal stress failure occurs in ply 1 at 122.5 lb. The ply 1 longi- 
tudinal (fiber direction) stresses are shown in figures 13 and 14, before and 
after dynamic equilibrium is reached by CODSTRAN. In figure 13, ply l longitu- 
dinal stresses at the center node are seen to reach their limit. In figure 14 
p y 1 has failed under the load and the corresponding longitudinal ply stress 
at the center of the plate is zero. The corresponding finite element mesh and 
displacement contours are shown in figure 15. For this example, ply failures 
that are initiated by the tensile stress limits usually activate the modified 
distortion energy failure criterion simultaneously. 


the plate in motion, 
is relatively slow, 
the end of the first 


. »f xt ’ the sajlie loading at the same rates is applied to the composite plate 
with the second type support where all boundary nodes are restrained from 
translation in all directions. Figure 16 shows the load-deflection relations 
for the three oading rates applied to the plate with the restrained supports, 
for the 0.1 Ib/psec loading rate, membrane action becomes much more important 
compared to the case of released supports and the displacements are signifi- 
cantly reduced as the restrained boundaries do not allow the majority of the 
plateto participate in transverse motion. Soon after the initial setting of 

inertial forces become negligible because the loading rate 
CODSTRAN analysis shows a 0.07 in. center deflection at 
load increment, under 30 lb load. During the next load 
increment, damage is initiated at the center node in ply l by reaching the 
stress limit in the transverse tensile stress. The 60 lb load that initiates 
damage is the same as what was required to initiate damage under the released 
type support condition. However, the displacement under the load is only two- 
thirds of the displacement under the released support case. Figure 17 shows 
t e transverse stresses that cause damage at the center node in ply 1. Fig- 
ure 18 shows the same stresses at much reduced levels immediately after damage 
when dynamic equilibrium is reached with the new structural properties that 
take the damage into account. Figure 18 indicates zero transverse stress in 
ply 1 at the center node where the damage occurred. Figure 19 shows the 
e ormed finite element model with values of the displacements after the ini- 
tial damage. Damage growth continues gradually as the transverse stress limits 
are exceeded in the other plies. Ply 2 fails at 180 lb, ply 3 fails at 210 lb 
and Ply 4 fails at 240 lb. All ply failures at this stage are in the trans- 
verse directions. Figure 20 shows the deformed plate with displacement con- 
tours after all plies are damaged in their transverse directions. 


The restrained boundary conditions affect the structural response at 
higher loading rates as well. With the restrained supports, initial damage 
occurs at lower magnitudes of the applied load since stresses are built up more 
quickly. However, the qualitative effect of the loading rate is similar in 

that the higher the rate of loading, the larger is the applied load to initiate 
damage . 


Figure 21 shows the deformed finite element model immediately after the 
initiation of damage, under 87.5 lb at 3.5 lb/psec. Figure 21 is comparable to 
figure 10 which shows the deformed structure at the time of initial damage by 
the same load under the released support conditions. Failure modes are also 
comparable as the first, second, and third plies fail simultaneously in their 
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transverse directions. However, there are some differences. First, the initi- 
ation of damage under the restrained boundary conditions starts at a lower 
load. Second, the damage region is spread over a much larger area of the lami- 
nate under the restrained boundaries. Figure 22 shows the transverse stress 
contours in ply 1 immediately before equilibrium iteration by CODSTRAN at the 
time of damage initiation. Approximately one third of the entire plate is par- 
ticipating in damage. Whereas in the case of the released supports, for which 
initial failure stresses are depicted in figure 11, initial damage was limited 
to the center node. For the restrained support condition, figure 23 shows much 
lower transverse stresses in ply 1 after CODSTRAN has progressed fracture until 
the structure reaches dynamic equilibrium. 

Additional loading after initial damage causes the fourth ply to fail in 
the transverse direction and the first ply to fail in the longitudinal direc- 
tion at 105 lb. 

Table I summarizes the initial damage and damage growth modes for all six 
loading cases considered. 


CONCLUSIONS 

An example structure is used to demonstrate the simulation of composite 
behavior under different loading rates, with two different boundary support 
conditions. At a relatively slow loading rate the dynamic response behavior 
is similar to that of a static loading; as the failure mode under slow dynamic 
loading has all characteristics of the failure that would also occur under a 
statically applied loading. This static failure mode is characterized by large 
deflections involving most of the structure. When the loading rate is suffici- 
ently high, inertial forces and local effects become significant and the fail- 
ure mode changes to a dynamic failure mode. In the dynamic failure mode 
inertial forces are significant, structural deflections are small, and large 
deformations are localized to the impact zone. At high loading rates, overall 
structural deformations at the time of damage initiation do not appear to be 
affected by the support conditions at the boundaries. However, stresses and 
damage propagation are influenced significantly by the boundary conditions. If 
the supports are restrained, damaged zones are much more widespread. On the 
other hand if supports are not restrained damage zones are localized to a small 
region under the applied load. 

With the presented research, a new dynamic computational capability has 
been implemented in CODSTRAN. The code is now capable of simulating the behav- 
ior of composites under any loading, geometry, material combinations, fiber 
orientations, tailoring, and boundary conditions. Structural degradation, 
delamination, fracture, and damage propagation are included in the simulation. 
The computational model is continually updated during the simulation of compos- 
ite response, using the degraded properties as local failures occur. Under 
impact loading, CODSTRAN can be used for the prediction of the type of failure 
as either localized impact penetration or more distributed flexural failure. 
CODSTRAN is also capable of monitoring degradation of structural properties 
such as natural frequencies, vibration mode shapes, and buckling loads during 
progressive fracture (ref. 3). Current research is aimed to verify and further 
develop the computational simulation by the addition of a new module for com- 
puting the global fracture toughness of composite laminates. Objectives for 
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continuing research include combining structural shell elements with solid ele- 
ments for more refined representation of the loading zone, computation of the 
impact loading by including the impacting object in the solution, representa- 
tion of strain rate effects in the constitutive relationships, fiuid composite 
structure interaction, dynamic impact loading and hydrodynamic shock loading 
of submerged composites composite degradation under cyclic loading, and non- 
destructive evaluation of remaining service life in composites. 
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END DISPLACEMENT 

FIGURE 1. - LOAD-DISPLACEMENT RELATIONSHIP UNDER STA- 
TIC AXIAL LOADING. 



FIGURE 2. - RECTANGULAR PLATE MODEL WITH NORMAL LOAD. 
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FIGURE 7. - STRESS CONTOURS AT DYNAMIC EQUILIBRIUM AFTER INITIAL 
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LIGURE 8. - DEFORMED PLATE AFTER INITIAL DAMAGE. 


LOADING RATE - 0.1 LBS''MICROSEC . 
APPLIED LOAD - 240 LBS. (A.E.) 


DISPLACEMENTS (IN.) 


.654- A 



FIGURE 9. - DEFORMED PLATE AFTER DAMAGE GROWTH. 
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FIGURE 10. - DEFORMED PLATE AFTER INITIAL DAMAGE. 
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FIGURE 11. - TRANSVERSE STRESSES AT INITIAL DAMAGE. 
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FIGURE 12. - TRANSVERSE STRESSES AT DYNAMIC EQUILIBRIUM AFTER 
INITIAL DAMAGE. 
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FIGURE 13. - LONGITUDINAL STRESSES AT DAMAGE GROWTH. 
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FIGURE 14. - LONGITUDINAL STRESSES AT DYNAMIC EQUILIBRIUM AFTER 
DAMAGE GROWTH. 
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FIGURE 15. - DEFORMED PLATE AFTER DAMAGE GROWTH. 


FIGURE 16. - LOAD-DISPLACEMENT RELATIONSHIPS FOR RE- 
STRAINED SUPPORTS. 
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FIGURE 17. - TRANSVERSE STRESSES AT INITIAL DAMAGE. 


LOADING RATE - 0.1 LBS^MICROSEC . STRESS (PSI) 

APPLIED LOAD - 60 LBS. (A . E . ) 



TRANSVERSE STRESSES FOR PLY 1 ^ 

FIGURE 18. - TRANSVERSE STRESSES AT DYNAMIC EQUILIBRIUM AFTER 
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FIGURE 19. - DEFORMED PLATE AFTER INITIAL DAMAGE. 


FIGURE 20. - DEFORMED PLATE AFTER DAMAGE GROWTH. 
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FIGURE 21. - DEFORMED PLATE AFTER INITIAL DAMAGE. 
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FIGURE 22. - TRANSVERSE STRESSES AT INITIAL DAMAGE. 


LOADING RATE - 3.5 LBS'MICROSEC . 
APPLIED LOAD - 87.5 LBS. (A . E . ) 


STRESS (PS I) 



FIGURE 23. - TRANSVERSE STRESSES AT DYNAMIC EQUILIBRIUM AFTER 
INITIAL DAMAGE. 
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